Incomplete fusion versus breakup competition with weakly bound projectiles 
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The importance of the breakup channel in the vicinity of the Coulomb barrier (1 < Ec.m./Vb < 2) 
is investigated for the medium weight ^Li + ^®Co system. Three-body final-state analysis of light- 
particle coincident data was carried out to disentangle, for the first time, the breakup contributions 
from other competiting mechanisms, a — d angular correlations show incomplete fusion components 
as significant as that from breakup process. Their strong coupling to total fusion is discussed within 
a comparison with predictions of continuum discretized coupled-channel calculations. 

PACS numbers: 25.70Jj, 25.70Mn, 25.70.Gh, 24.10.Eq 



The study of fusion reactions in the vicinity of the 
Coulomb barrier provides a fascinating challenge for the- 
ories of quantum tunneling leading to an irreversible com- 
plete fusion (CF) of the interacting nuclei into the com- 
pound nucleus (CN) HQ. The fusion probability is sen- 
sitive to the internal structure of the interacting ions as 
well as to influence of the other competing mechanisms 
such as nucleon transfer and/or breakup (BU) which are 
known to affect the fusion. The fusion cross section en- 
hancement generally observed at sub-barrier energies is 
understood in terms of dynamical processes arising from 
couplings to collective inelastic excitations of the target 
and/or projectile. However, in the case of reactions where 
at least one of the colliding ions has a sufficient low bind- 
ing energy so that BU becomes an important process, 
conflicting experimental |^ ^ |^ and theoretical results 
are reported ^,J^ 8, 9]. 

A great experimental effort, involving both (loosely 
bound) stable and unstable nuclei, has been devoted to 
investigate the specific role of the BU channel 0, 0, . 
The weak binding of these systems can also lead to in- 
complete fusion (ICF)/ transfer (TR) processes playing 
an important role. Several attempts to clearly identify 
ICF in ^Li and ^Be induced fusion with fissile targets 
0, 0| and medium-mass targets 0, 0| have been 
made. CF requires the formation of CN containing all 
the nucleons of both projectile and target. If only part 
of the projectile fuses with the target, with remaining 
fragment emerging from the interaction region, then ICF 
is defined (in this case, the BU process is followed by 
fusion) A tranfer-reemission process may lead to a 
final state similar to ICF [l7j . 

The recent availabilit y of lig ht-mass radioactive ion 
beams such as ^He [H El US 113. "Be [H, and "F 
|23|. and the renewed interest on reactions involved in 
astrophysical processes motivated the investigation 
of fusion reactions involving very weakly bound and/or 
halo projectiles around and below the Coulomb barrier. 



Clearly a full understanding of the BU process and its ef- 
fects on near-barrier fusion is fundamental in order to be 
able to understand the dynamics of reactions involving 
radioative nuclei. This requires systematic and exclusive 
measurements covering a wide range of processes, sys- 
tems and energies. We choose to study both the total 
fusion 25] and BU ^26] of ^'^Li with the intermediate- 
mass target ^^Co. 

In this Letter we address the competition between the 
several reaction processes and CF. A three-body kine- 
matics analysis, in which we separate the contribution of 
the BU and the ICF, is presented for the first time. 

The experiments were performed at the University of 
Sao Paulo Physics Institute. The ^Li beam was delivered 
by the 8UD Pelletron accelerator with energies Ejof, = 18, 
22, 26 and 30 MeV, and bombarded a 2.2 mg/cm^ thick 
^^Co target. The detection system consisted of a set 
of 11 triple telescopes [l^l separated by 10°, for which 
light particles can be detected with a very low-energy 
threshold (0.2 MeV for d and 0.4 MeV for a particles). 

In this work we will concentrate on a ~ d coinci- 
dences, which are usually fully attributed to a "BU 
process" 0. Depending on the angular combination, 
we observe well defined peaks, in the total kinetic 
energy spectrum, corresponding to the sequential BU 
of ^Li in its first excited state with E* = 2.19 MeV. 
No other discrete excited states are observed. When 
analysing the a + d coincidence yields from the ^Li 
induced reaction, we have to consider the contributions 
of other processes than BU, leading to the same parti- 
cles in the final state. The processes to be considered are: 



i) ^Li + 59Co - 

ii) 6Li + 59C0 - 

Hi) ^Li + 59Co 
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Process i) is identified as the sequential BU of ^Li. The 
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final state can also be reached through a direct BU. Pro- 
cess ii) can be identified as incomplete fusion of d + ^^Co 
{d ICF) with the subsequent reemission of a deuteron 
from the excited ^'^Ni. This process could also be consid- 
ered as a d transfer followed by a d reemission from the 
^^Ni nucleus. The same observations are valid regarding 
process Hi), for which cither incomplete fusion of a -|- 
^^Co {a ICF) or a transfer could occur and an a particle 
is reemitted. Process iv) corresponds to the a — d sequen- 
tial decay of the ^^Zn CN. The contribution of CN decay 
is considered negligible, as confirmed by predictions from 
statistical model codes. 

In order to investigate the competition of the above 
processes, we performed a complete three-body kinemat- 
ics analysis. As we know the detection angle and energies 
of both the d and the a particle, from the three-body 
kinematics equations we can determine the energy and 
emission angle of the remaining ^^Co nucleus and from 
this, all the quantities of interest, such as Q- values and 
relative energies. By generating Q-value spectra, we ob- 
serve from the data, for most of the events in the final 
state a + d + ^^Co, the residual "''''Co nucleus mostly in 
its ground state, in which the events were gated. Prod- 
ucts from the sequential ^Li BU a + d, are focused inside 
an angular cone. If we assume that process i) is occur- 
ring, the relative energy E^-d between a + d and d in 
the rest frame of ^Li can be calculated. By fixing for 
instance, the detection angle of the a + d particle and 
varying the d detection angle, we can follow the behavior 
of the relative energy E^-d. In Fig. 1(a) is shown, for 
Fiiab = 29.6 MeV (corrected for the energy loss to the 
center of the target), the relative energy ^a-d as a func- 
tion of the deuteron detection angle Od, for a fixed angle 
6a = 45° of the a particle. It is interesting to notice that, 
within the angular range where the sequential BU of ^Li 
in its first excited state is kinematically allowed (delim- 
ited by the two vertical dashed lines) , the relative energy 
Ec_d is constant. The ^a-d value is consistent with the 
sequential BU of ^Li in its first excited state. Outside 
this region it is no longer constant, which suggests the 
presence of other processes. In addition to ICF or TR, 
a direct BU occuring close to target with strong nuclear 
field could also give values where the quantity E^-d is 
not constant. 

In Fig. 1(b) we show the behavior of the laboratory 

a particle kinetic energy Eq, as a function of 9^ for a 
fixed angle = 45°, and for Eiab = 29.6 MeV. Here, E„ 
is taken as the centroid of the experimental coincidence 
spectra. We observe that the average energy Eq, is con- 
stant (independent of the momentum of the deuteron), 
except in the angular range where the sequential BU of 
®Li is present. Considering now process ii), if this bi- 
nary process is occurring with an intermediate stage, for 
a given angular combination 9a and 9d the energy E„ 
is uniquely determined once the ^'^Ni excitation energy 
is defined. We can then conclude from Fig. 1(b) that 
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FIG. 1: (a) Relative energy E^-d as a function of the detec- 
tion angle 6d- (b) Average a energy Eq, as a function of the 
deuteron detection angle 8d- (c) d energy in the rest frame of 
the decaying ^'^Ni as a function of 6d- The angles with two 
experimental points correspond to two kinematical solutions 
for the sequential BU. 



if process ii) is dominant over the angular range where 
Ea is constant, the average ^^Ni excitation energy is also 
constant, and in this case Eei^^ = 25 MeV. This is a very 
important result as a clear indication of a CN-type reac- 
tion. This constant value is consistent with the assump- 
tion of an ICF d + ^^Co where the d has the projectile 
velocity. As a consistency check, we also calculated the d 
energy in the rest frame of the decaying ^^Ni. If the ^-"^Ni 
is in equilibrium and the excitation energy is constant 
over the mentioned angular interval, the d energy should 
be constant and independent of the d emission angle 9d- 
This is indeed observed, and is shown in Fig. 1(c). The 
same behavior is observed for other configurations with 
different fixed angles 9a- On the other hand, if a TR is 
assumed, an optimum Q-value can be calculated and in 
this case the residual nucleus ^^Ni would have an excita- 
tion energy E|i^. — 22 MeV. This value is also close to 
the experimental one. Taking into account the different 
relations available for the optimum Q-value calculation, 
this could also be closer to 25 MeV. 

The process defined in Hi) could also be present in the 
coincidence data. In order to check that, the same analy- 
sis described above was performed for the situations with 
a fixed angle 9d as a function of 9a- A similar behavior is 
obtained, and this leads us to conclude that, within our 
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FIG. 2: (Color online) Experimental in-plane angular corre- 
lations and the predictions for the contribution of d ICF and 
a ICF. The contribution of the BU is also shown. The total 
sum is shown as a solid line. 



sensitivity, there is also a significant contribution from 
process Hi). 

In order to quantify the contribution from processes 
j), m) and Hi), we constructed the a — d angular corre- 
lation functions through calculation of the double differ- 
ential cross sections a / {dfladQd) ■ The absolute cross 
sections, the product of the number of particles in the 
target per unit area and number of particles in the beam 
{NaNb) for each run, was calculated and normahzed to 
the elastic scattering data we measured. The uncertain- 
ties in the experimental points (about 10% to 40%) are 
due to statistics, the determination of NaNb and, the 
geometrical determination of the detector solid angles. 

As mentioned above, a priori, events from TR and ICF 
process are indistinguishable due to the fact that the in- 
termediate nucleus, in both cases, is populated in the 
continuum, and a statistical description for these pro- 
cesses seems adequate 28,j. The model we propose as 
follows is based on an ICF picture. 

Assuming that for the d ICF (a ICF) process, the ex- 
cited ^^Ni (^^Cu) is in equilibrium prior to the d (a) 
reemission, a model can be utilized to describe the de- 
cay. We choose here the approach developed by Halpern 
[29j. It consists of a classical model for emission from 
a spherical rotating nucleus. If a rotating nucleus is as- 
sumed to rotate around an axis that is perpendicular to 
the reaction plane, then the angular distribution of the 
evaporated particles is isotropic in the equatorial plane. 
Due to the centrifugal force, the yield is concentrated in 
the equatorial plane and decreases toward the poles. If 
the rotating axis is normal to the reaction plane, then 
the equatorial plane is the reaction plane. 

The yield of evaporation particles as a function of the 



polar angle tp defined with respect to the rotation axis is 
given by: 



(1) 



where Yq is a normalization factor and X is the ratio of 
rotational kinetic energy to the thermal nuclear energy: 



X = 0.5(J+ i)V2/T 



(2) 



where J is the spin of the rotating nucleus, / — ^R^ is 
the moment of inertia, and the temperature T can be 
estimated from E* = aT^, with a being the level density 
parameter. 

The angle between the rotational axis and the z axis 
which is normal to the reaction plane, is defined to be 7. 
The angle 7 is assumed to be gaussian distributed. From 
Eq. (1), the angular distribution is: 



W{9,^)^ / d7 cxp(-7V27o)>"W 



(3) 



with cosijj = COS7COS0 -I- sin7sin(^sin0. 

The angular distribution W{6, (p) is the same as the 
experimental quantity d^a/{diladfld) in the rest frame of 
the rotating nucleus. In the laboratory reference frame, 
W{9, (j)) is centered at the recoil direction in the primary 
process. It is important to remark that by using this 
model, the angular correlation cj) dependence is also taken 
into account. Therefore, the integration is performed in 
a better way than assuming an isotropic (f> dependence. 

For d and a emission, the best values for Yq, X and 
7o are obtained from fits to the angular distributions 
provided by the statistical code STATIS [Ull^l. For 
the calculations, a fusion process d + ^^Co — > ^-'^Ni* or 
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Co 
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Cu* was assumed, with a bombarding 



energy forming the excited CN in an excitation energy 
corresponding to the most probable value observed ex- 
perimentally. 

In Fig. 2(a) we show the a — d angular correlation for 
a fixed 9a = 45° and Eiab = 29.6 MeV, together with the 
model predictions for the d ICF and a ICF. The shape 
of the a ICF correlation is obtained from the model pre- 
dictions for the angular correlations with 9d fixed. In the 
same manner, the d ICF shape presented in Fig. 2(b) 
is obtained from the predictions for the angular correla- 
tions with 9a fixed. The sequential and direct BU con- 
tribution is also shown in Fig. 2. It is obtained through 
the subtraction of the incoherent sum of the d and a 
ICF contributions, from what would be the best fit to 
the data. The two peaks lying around 45° correspond to 
the sequential BU of the 2.19 MeV unbound state of ^Li. 

By numerically integrating the angular correlation 
function in 9 and 0, for each process, the differential cross 
section da/dQa (dcr/dOd) is obtained for the fixed 9a 
{9d)- The procedure is repeated for all the other angular 
correlations with different fixed 9a or 9d. 
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analysis of a — d coincidences of ^Li + ^^Co at several 
near-barrier energies, we observed a significant contribu- 
tion from the ICF process, with a cross section compa- 
rable to the one from BU process. This suggests that 
ICF should be accounted for in the coupled channels cal- 
culations to explain the fusion of weakly bound nuclei 
inhibition and/or enhancement at near-barrier energies. 

The authors thank FAFESF and CNPq for the finan- 
cial support and also Frof. C. Signorini and Frof. B. R. 
Fulton for the careful reading of the manuscript. 



FIG. 3: (Color online) ®Li -I- ^^Co excitation functions for d 
singles, d ICF, a ICF, BU and fusion :25|. The BU CDCC 
prediction is shown as a solid line. The total fusion CDCC 
prediction is shown as a dashed line. 



The BU cross section is obtained with the utilization 
of the d singles cross sections. We assume that the d 
singles spectra present predominantly the contribution 
of BU, d ICF and a ICF. The total BU cross section 
for a given bombarding energy is then obtained by the 
difference between the total d singles cross section and 
the total d ICF -I- a ICF cross sections. 

The exclusive BU cross sections for the resonant states 
in ^Li have been calculated by the CDCC formalism 
using a cluster folding model with potentials that de- 
scribe well the measured elastic scattering angular distri- 
butions [2g. The CDCC calculations for °Li were per- 
formed with the code FRESCO assuming an a + d cluster 
structure, similar to that described in |3l| . The binding 
potentials between a + d were taken from [s^ and the 
a + d continuum was discretized into series of momen- 
tum bins of width Sk = 0.2 fm"-"^ (up to fc = 0.8 fm"-'^) 
for L = 0,1,2 for ^Li, where Tik denotes the momentum 
of the a + d relative motion. All couplings, including 
continuum-continuum couplings, up to multipolarity A 
= 2 were incorporated [S^l- The total calculated BU 
cross-sections for ^Li were obtained by integrating con- 
tributions from states in the continuum up to 11 MeV. 
The results of full CDCC calculations are displayed in 
Fig. 3. The solid line corresponds to the final results of 
the BU CDCC calculations, and the symbols represent 
the experimental excitation functions obtained for the d 
singles, and for the d ICF, a ICF and the BU processes. 
The d ICF and a ICF cross sections are comparable to 
the BU cross sections for all the bombarding energies. 

In summary, we have presented a three-body kinemat- 
ics analysis of coincidence measurements with the aim of 
disentangling the contribution of BU and competing pro- 
cesses in reactions with weakly bound nuclei. From the 
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